Noise in an SSET-resonator driven by an external field 
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We investigate the noise properties of a superconducting single electron transistor (SSET) cou- 
pled to an harmonically driven resonator. Using a Langevin equation approach, we calculate the 
frequency spectrum of the SSET charge and calculate its effect on the resonator field. We find 
that the heights of the peaks in the frequency spectra depend sensitively on the amplitude of the 
resonator oscillation and hence suggest that the heights of these peaks could act as a sensitive signal 
for detecting the small changes in the amplitude of the drive. The previously known results for the 
effective amplitude-dependent damping and temperature provided by the SSET for the case of a 
low frequency resonator are generalized for all resonator frequencies. 

PACS numbers: 85.85.-l-j, 85.35.Gv, 74.78.Na,73.23.Hk 



CD 



a 



o 



> 

^^ 

\>- 

o 

00 

o 



X 



Over the past few years experiments on superconduct- 
ing circuits have produced some rather impressive results. 
Superconducting elements have been used to produce 
two-level systems of various kinds which can be consid- 
ered as artificial atoms [E S S S] j and superconducting 
stripline resonators can act analogously to optical cav- 
ities for microwave fields Q. This field of study has 
been referred to as circuit QED in analogy with cav- 
ity QED. One advantage of superconducting circuits is 
that rather than simply investigating the behavior of the 
system through the field emitted or reflected by the cav- 
ity, other elements can provide additional information. 
Mesoscopic conductors coupled to the resonator can pro- 
vide additional information and at the same time the 
back action may lead to some interesting and subtle dy- 
namics. Such back-action dynamics have also received 
considerable attention in the context of a mesoscopic con- 
ductor used to investigate the behavior of a mechanical 
resonator [1, 0, S i, ES [O, El [11 • 

One system that is of particular interest is that of a su- 
perconducting single electron transistor (SSET) coupled 
to a resonator, either mechanical [Tj, [l^, ]1Q, \VA\l&l . UM . 
^0] or composed of a superconducting stripline [23] . The 
coherent transport through this device at the Joscphson 
quasiparticle resonance [23jl23 (JQP) allows a very low- 
noise current, and at the same time the sensitivity of the 
SSET to charge means that the resonator-SSET coupling 
is significant. The SSET biased at the JQP resonance can 
be considered analogous to a three-level atom, and the 
coupled system therefore shows behavior related to that 
of a micromaser |2l|, [23 • 

The low noise and relatively strong coupling in this 
device allows the observation of the non-trivial coupled 
dynamics that arises from the interaction between the 
resonator and SSET. When the SSET is biased so as to 
absorb energy from the resonator, the SSET acts as an 
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additional source of damping, and the resonator can be 
cooled below its physical temperature. This modification 
of the steady state of the resonator due to a back-action 
induced effective damping and temperature has been ob- 
served for an SSET coupled to both mechanical and su- 
perconducting resonators [l9|, |23|. When the SSET is 
biased to the other side of the JQP resonance, energy 
is transferred from the SSET to the resonator, leading 
to an effective negative damping. For a superconducting 
stripline resonator, the quality factor can be high enough 
that the negative damping dominates and the resonator 
is driven into an oscillatin g st ate even in the absence of 
an external drive |14l . [isl . [2q . This has recently been 
observed for a superconducting stripline resonator cou- 
pled to an SSET, and the field transmitted through the 
resonator measured [2^. When the resonator is in this 
oscillating state, the SSET has a periodic but not har- 
monic response, with the details of its dynamics depend- 
ing rather sensitively on the amplitude of the resonator 
oscillation. In order to probe this behavior, we must go 
beyond the time-averaged mean field behaviour [l5| and 
investigate the correlations present in the system. 

With a stripline resonator, as opposed to a mechanical 
one, there are two independent probes of the system - 
the field in the superconducting stripline, and the cur- 
rent that passes through the SSET. This then allows the 
investigation of the system by two different experimen- 
tal methods. The noise properties of the SSET-cavity 
system are transferred to both the current and reflected 
field, so we can consider these to be probing the noise 
present in our system. The effect of the thermal noise in 
both the SSET and the cavity can be made very small, 
so the noise detected largely arises from the finite level 
structure of the SSET, and is in that sense is a purely 
quantum noise. 

In this paper, we consider an SSET that is coupled to 
a resonator that is driven into an oscillating state. We 
calculate the charge noise on the SSET and the effect this 
has on the cavity field. We also consider the case when 
the resonator's motion depends on the effective damping 



and noise arising from the SSET, using a linear response 
approach ^] . We go on to show how the response of the 
system at the sidebands could be used to detect small 
changes in the driving amplitude. Although we focus our 
analysis on the case of a superconducting stripline, which 
has a higher oscillation frequency, our main results are 
also valid for the lower frequency mechanical resonators. 
In our calculation, we take advantage of the fact that the 
resonator oscillations change on a timescale that is slow 
compared to the SSET timescales, which is valid as long 
as the resonator damping and coupling to the SSET are 
weak. This means that the SSET response is essentially 
the response to a purely harmonic drive. This approach 
has previously been used for both mechanical and op- 
tomechanical systems [1^ [2^|. However, by including a 
fluctuating Langevin term, we can go beyond the previ- 
ous mean field results to calculate the noise. The effect 
of this response on the resonator can then be calculated. 

In Sec. U we introduce our model of the system and 
Langevin equations for the resonator field and the SSET 
variables. In Sec. |TT] we introduce the main approxi- 
mations made and solve the Langevin equations for the 
SSET under the influence of the periodic driving pro- 
vided by the resonator. Section IIIII describes the fre- 
quency spectrum of the charge on the SSET and gives 
approximate expressions valid when the resonator fre- 
quency is large. In Sec. IIVI we show how the charge 
dynamics affect the field in the cavity. We derive expres- 
sions for the effective SSET damping and temperature of 
the resonator, going beyond previously published results 
to derive expressions valid for all resonator frequencies. 
Section |V] describes how the response of the system at 
multiples of the resonator frequency could be used to de- 
tect small changes in the driving amplitude and Section 
IVII describes how the field in the cavity could be mea- 
sured through a transmission line coupled to the cavity. 
In Sec. IVIII we present our conclusions and in appendix 1^ 
we present more details of the derivation of the Langevin 
equations for the SSET. 



MASTER EQUATION AND LANGEVIN 
EQUATIONS 



The master equation for a superconducting single elec- 
tron transistor biased at the Josephson quasiparticle res- 
onance and capacitively coupled to a resonator [Ij, [l^ , 



P — Cp ~ — — [Hco, P] + CleadsP + ^dampingP- (1) 



The effective Hamiltonian is given by, 

i/co = -?iAcr22 - hej (cro2 + Cr2o) 



+ a^e-'"^°*) 



— huJn 



2x, 



■(a + at)((7ii + 2CT22) 



(2) 



where the operators aij = \i){j\ represent operators on 
the SSET island charge states |i), A is the energy differ- 
ence between states |2) and |0) and ej is the Josephson 
energy of the superconductor. The frequency of the res- 
onator is ujo and Ajj , lud give the strength and frequency 
of the external driving respectively. The resonator-SSET 
coupling is described by the parameter Xg which mea- 
sures the shift in the equilibrium field of the resonator 
brought about by adding a single electronic charge to the 
SSET island [iil, and cc, = (?i/2eoWo)^- The tunneling 
of quasiparticles from the island is described by 



C-leadsP = --^ [{0'22 + Crii,/9}_|_ 

-2 ((712 + CToi) P (o-21 + CTio)] 



(3) 



where we have neglected the position dependence of the 
tunnel rates as being of lesser importance than the co- 
herent coupling |14l . llSJ I . This simplification means that 
the dissipation takes a Lindblad form, and also means 
that the master equation is essentially equivalent to that 
of a resonator coupled to a double quantum dot. The 
dissipation and fluctuations arising from the resonator's 
surroundin gs a re described by the usual quantum optical 
expression [S^l, 

CdampingP = - -^{n + I) {a^ ap + pa^ a - 2apa^) {4) 
—n (aa) p -f paa) — 2a* pa) 

where this form for the dissipation is valid as long as 
the dynamics of the system are slow compared to the 
correlation time of the bath. 

An alternative description of the system is given by a 
set of Langevin equations describing the coupled system. 
This gives an equivalent description of the first and sec- 
ond moments of the system, which is what is required for 
a noise calculation. We shall see later that this Langevin 
form is convenient for dealing with the external drive, as 
a transformation can be made which allows us to solve 
the equation analytically. 

The Langevin equations are equal to the semiclassical 
equations plus a fluctuating noise term. For the resonator 
operator a we have. 



-ibj^a ■ 



le 



-a - a^e-*""* + T^a 



consists of a coherent part, described by the Hamilto- 
nian, Hco, together with two dissipative terms Cieads and 
^damping which dcscribc quasiparticle decay from the is- 
land and the surroundings of the resonator respectively. 



+iu!o—^{aii -1-20-22) 

ZXf, 



(5) 



where r]a represents the standard white noise term on 
the resonator defined by (?7a(i)) = , {'i]i{t)Va{t')) = 



<5(i-i')7e.n and {va{tH{t')) = S{t-t')j,,{n + l). This 
couples to the equations for the SSET charge operators 
^u = l*)OI UM- The Langevin equations for the projec- 
tion operators are given by, 

<^oo{t) = iej{ao2{t) - (T2o{t)) +Taii{t) +r]oQ{t). (6) 
<yii{t) = ra22{t)-Tan{t)+T^iiit) (7) 

^22 W = -iej{ao2{t)-<72oit))-Ta22{t)+V22{t)i8) 

which also couple to the equation of motion for the off- 
diagonal terms describing the coherence between levels 
|0) and |2), 



O'02(i) 



-«ej(o-22(i) - cro(){t))--ao2{t) 



(9) 



-i{A^ 



UJqXs 



{a{t)+a\t)))ao2{t)+V02{t). 



We need to find the properties of the correlators for the 
SSET noise operators rjij. The Langevin equations, in- 
cluding the properties of the noise operators, can be eas- 
ily derived from the master equation by requiring that 
the two forms give the same equations of motion for the 
second moments of the operators, determining the corre- 
lation properties of rjij (see appendix [X)) . For example, 
for the operator representing the noise on cto2 we have, 

(^02(i)r/2o(i')> = r((aoo(i)) + (aii(t)))5(i-i')-(10) 

There are a few points to note about the noise operators 
for the SSET. Firstly, in the limit ksT < hej, the SSET 
experiences no direct thermal noise. The "noise" terms 
7702 (t) therefore arise solely from the finite level structure 
of the SSET and the resulting commutation relations, 
and in that sense can be considered to be purely quantum 
noise. Secondly, we are giving an approximate treatment 
of this quantum noise by requiring that the first and sec- 
ond order moments calculated by the quantum Langevin 
equations are equivalent to the first and second order mo- 
ments as determined by the master equation. In general 
these would not suffice to determine all the moments in 
the system. Finally, we note is that we are describing a 
system under periodic driving which will therefore tend 
to a periodic behavior in the long time limit, rather than 
a fixed point. In a finite level system, the correlators in- 
volve the average values of the operators such as (crii(t)) 
which depend on t in a driven system. This means that 
we have time-dependent noise correlators for the SSET 
noise operators. 

The equations for the operators cri2 ,crio decouple from 
the others, so Eqs. ([Ml]), along with their correla- 
tors, completely determine the dynamics of the SSET- 
resonator system, up to second order. 



II. SOLVING THE ao2 LANGEVIN EQUATION 

In this section we solve the Langevin equation (Eq. 
([5])) for the off-diagonal charge operator cro2- In order to 



progress we make two assumptions: first that the Joseph- 
son energy is rather small, and second that the ampli- 
tude of the resonator changes slowly compared to the 
incoherent dynamics of (T02, i-e. that the total damping 
due to both the resonator environment and the SSET 
It = 7ex + 7SS is much smaller than the quasiparticle 
decay rate F. 

In the limit that the Josephson energy is much weaker 
than the quasiparticle decay, ej ^ F, the occupation of 
the charge states cth, CT22 ^ 1, and the equation of mo- 
tion for (T02 decouples from the other charge equations. 
We are assuming the resonator amplitude can be treated 
as constant on timescales relevant to the SSET dynam- 
ics, an approximation which has proved useful in related 
mechanical and optomechanical systems [1^, [2^|. Our 
derivation closely follows these methods, but also incor- 
porates the fluctuations described by the 77 terms. We 
replace the term describing the field in Eq. ([9]) with a 
cosine oscillation at the driving frequency with magni- 
tude A i.e. (a(i) -|- a^t)) « Acos(w_Di). For the case 
of weak back-action damping, the amplitude is simply 
given hy A = Ad/Hcoo — ojdY + llxY^'^- When the back 
action damping is significant, the amplitude must be de- 
termined self-consistently - an oscillation of amplitude A 
will be stable if the resulting total damping is zero. 

This then means that the field simply appears as an 
harmonic drive acting on the SSET, and Eq. ^ be- 
comes, 



(T02(t)= iej- 



F uji)Xs \ , \ / \ 

--zA-i AcOSUJut <T02(t)+T]02(t), 

2 x^ J 

(11) 



where the first term on the right hand side multiplies an 
implicit identity operator. 

If we make a transformation to eliminate the driv- 



ing term, CT02 = o'02e' 



-iz sin ujot 



, we can then find the 



Fourier transform of the transformed operator ao2 in 
terms of Bessel functions of the first kind, J„(z), where 

Z = {uJoXsA)/{uJDXq), 



1702 (w_f) 



i'^J Z) K'^F - UJDn)Jn {-z) 
n 

r/2 + i{ujF-A) 

^ J g-»^ptg-»2sina;Dt^g2(i)dt 



F/2 + i(wF- A) 



(12) 



and we see that 0-02 consists of a systematic response to 
the drive at multiples of the driving frequency wd, plus 
a noise term. 

Equation P^ gives the Fourier component in the 
transformed picture, so we convert back to the untrans- 



formed picture to obtain, 

iejS{ujF - iODn')Jn'-n i-z) J„ (z) 



(J02{uJf) = ^ 



T/2 + i{ujD{n' ~n)- A) 

OO 



2tt{T/2 + i{uoF - ujD-n - A)) 



(13) 



which gives an expression for the Fourier transform of 
f 02 (i) consisting of a mean field and noise- induced term. 
From this expression we can go on to calculate the spec- 
trum of the charge on the SSET and hence of the cavity 
field. 



III. NOISE ON THE SSET 

The frequency spectrum of two fluctuating terms 
f{t),g{t) is defined by, 



Sfg{LO,t) = J dre--(/(i + r)5(i)) 

— OO 

OO 

^ 27T [ dujFe'^^^+^')'{f{uj)g{uJF)) (14) 



where in the second line we have written the spectrum 
in terms of the Fourier transforms of the individual func- 
tions. Note that as defined, the spectrum includes corre- 
lations due to the systematic motion of the terms as well 
as due to the fluctuations. If the functions f,g are peri- 
odic then the above expression retains a periodic time de- 
pendence. We therefore define a spectrum averaged over 

■it/ujo 

a single driving period, Sfg{uj) = ^ / dtSfg{u>,t), 

—■k/uiu 

which is what we would expect to observe in experiment. 
Inserting Eq. (jl3|l into this expression and noting that 
in the small-ej limit we can approximate the time de- 
pendent correlator given in Eq. p^ as {r]o2{t)ri2o{t')) ~ 
r6{t-t'), gives, 



Oa 



20 (^) 



y 2iT5{Ld — uoD'n 

n 

+E 



E 



ejJn' i-z) Jn-n' [z) 



5 + i{u)Dn' - A) 



rjn (zY 



{uj-ujDn-A)^ + r'^/A' 



(15) 



The first term corresponds to the correlations arising 
from the systematic response of the resonator to the driv- 
ing force. For a purely monochromatic drive, the peaks 
are (5— functions, but in reality the (5— functions would be 
replaced by peaks with total power unity and a width 
determined by the linewidth of the driving. The second 



term describes the additional noise due to the finite level 
structure. 

We have calculated the noise on the off-diagonal ele- 
ment of the SSET density matrix, cto2- However, this 
does not enter directly into the current or cavity noise. 
Instead, we have to calculate the frequency spectrum of 
the total charge on the SSET, ace = cii + 2(722. 



A. Systematic Charge Oscillations 

Fourier transforming Eqs. ^S^ gives an expression 
for the charge in terms of (Jq2{ojf)- Inserting Eq. p^ 
into this gives an expression for the charge consisting 
of a systematic part, and a fluctuating part containing 
the noise operators ryo2,'7ii and 7722- We now go on to 
calculate the fluctuation spectrum of the charge due to 
the systematic and noisy terms. 

The systematic component of the charge, <Jq2{^f) is 
given by. 



<^!ci^F) 



(3r + 2iLUF)£] 
(F + iojfY ^, 



- '^^ Jri' {.z) S{ujf - i^on) 

^n'+n \Z } 'Jn'-n \Z } 



^-i{ujDn' + A) 



^ + i{ujDn' + A) 



(16) 



and we see that the mean field response of the charge is 
a series of (5— peaks at multiples of lod. 

The systematic motion of the charge will show up in 
the time correlation function and hence the fluctuation 
spectrum. We insert Eq. P^ into Eq. (|14p to obtain the 
part of the charge spectrum due to systematic evolution. 



siiuj) 



E 



E 



2'i:5{u} + nuj£i)e J 

Jn'(z) Jn+n' (z) 

E - z(wDn'-t-A) 



4 9r2 + Auj^n^ 



{T^+ujln^y 



(17) 






For F, A ^ ujf), this expression is dominated by the n' = 
term in the sum, and the systematic noise reduces to. 



Scci^) ~ V27r(5(w + ncjD)47 T^Jo{zf 



Jn {z) , J-n {z) 



-iA 



? + .A 



(18) 



The symmetry properties of the Bessel functions then 
mean that the bessel functions will either add or cancel 
depending on whether n is odd or even, so separating out 



the two cases we find, 



Siico) 



4ej Jo {zY ■^^ 2TrS{LLj + nujo] 



■E 



+ 4A2 



Jn {z) 



.2(l + (-l)") , .^2(l-(-l)") 



(19) 



B. Charge Spectrum Due to Fluctuations 

As well as the systematic response to the driving 
force, the charge spectrum also contains a part due 
to fluctuations. The Fourier transform of the charge 
c^cc = cTii + 2(722 contains the fluctuating terms corre- 
sponding to both the diagonal and off-diagonal opera- 
tors, 7711, 7722, '702- Inserting these into Eq. (|14p . we write 
the fluctuation- induced part of the charge spectrum as, 



nri _ Qdiag . Qri02,V20 i across 



(20) 



where S*^*"^ corresponds to the part arising from fluctu- 
ations on the diagonal operators 7711,7722, the term 5°^ 
arises from 7^02 , and S^™'"* from the correlations between 
these. Note that while {rio2{t)ri2o{t')) can be approxi- 
mated to a constant in the small-ej limit, other correla- 
tors give, for example, {r]22{t)r]22it')) = T{a22{t))S{t-t'), 
i.e. we have a time dependent function multiplied by a 
delta function. As (1722 (i)) is simply the systematic part 
of CT22, we can insert the Fourier transform of this into 
our calculation. After some algebra, we find. 






E 



JnizTTe] 



-t- (w_d71 -I- A)2 
JnizfTe^j 



(21) 



SF^ + 2cj2 

(F2+w2)2 

9r2 + 4c^2 

(F2+C.2)2Z_.I-+(^„^^„_A)2 

(22) 



The expression for the cross terms is rather unwieldy, but 
can be written as. 




i2r'^ + 4tj2 



2iLur 



(F2+C^2)2 



(23) 



Jnizfrej 



(^ - i(w_Dn- A))(^ + i{uj -ujon -A)) 



In Fig. [T] we plot the charge noise as a function of w, 
and see it consists of a series of peaks (of width ~ F) 
at integer multiples of the driving frequency on top of a 
background peak centered omo — Q. It is these sideband 
peaks that give information about the driving frequency 
and amplitude, and so we would like to find some simple 
expressions that tell us if these peaks can be observed. 
The peaks are more pronounced in the limit that the 
resonator frequency is larger than the quasiparticle decay 
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FIG. 1: (color online). Charge noise as a function of u for dif- 
ferent driving frequencies lod = (a)0.5, (6)2, (c)10 relative to 
the quasiparticle decay F. The different curves show different 
driving amplitudes 2 = 2,3, 4. The inset (d) in the u>d = 10 
plot shows that the large central peak is highly suppressed 
at Jo{z) — 0, here at z — 2.4048. The other parameters are 
ej = 1/16, r = 1,A= -0.5, 



rate, and this is what we would expect to be the case 
experimentally for superconducting stripline resonators. 

The spectrum simplifies considerably in the large w^i 
limit. The diagonal term, Eq. ([2T|) . consists of a peak 
at a; = 0, multiplied by a sum that is independent of w. 
We can approximate the sum by the 77. = term in the 
F, A ^ lod limit to obtain. 

The other terms describe a similar peak at lu — D, 



multiplied by a series of peaks. We wish to find an ex- 
pression for the heights of these peaks. In the F, A <C lod 
limit, we find that the heights of the noise terms at inte- 
ger values of the driving frequency are, 






4J„(z)2 



r^^. 



2 r2 






{ujonf ri^^2 



(24) 



and we note that the cross terms are negligible in this 
limit and can be neglected. Thus we can write the height 
of the charge peaks as 



S'ilci^Dn) 



(4J„(z)2 -f 2Jo(z)2 



r^} 



2T:{LUDny 



^ + A2' 



(25) 



The visibility of the peaks depends on the contrast be- 
tween the noise at the peaks and the noise between the 
peaks [37| . 

Between the peaks, two terms in the sum are relevant 
which gives, to leading order. 



1. 



SliuD{n+-)) 



2Mzf 



^e'j 



(26) 



The ratio reaches a minimum when Jn{z) = 0, at which 
point we see the ratio is {n + i)^/n^. The ratio becomes 

large when Jo{z) = 0, and is of order uj^n? /{^ + A'^). 
Examples of the on- and off-peak heights as a function 
of the driving are shown in Fig. ^ , along with the ap- 
proximations to these given in Eqs. (|25II26[) . 
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IV. CAVITY FIELD 

We now consider the effect of the SSET on the cavity. 
We recall the Langevin equation for the cavity field. 



lex . . 

—iUJr,a a — lA 

2 

+iuJo—^{aii + 2(722) 



oe '^"'+Va 



(28) 



where A]j represents the amplitude of the classical driv- 
ing field. 

If we neglect the back-action as weak, then the steady 
state amplitude A = \a\ of the resonator is simply given 

by, 



A = 



A 



D 



{ujQ - uJoY + ^ 



(29) 



FIG. 2: (color online). Comparison of the heights of the peaks 
at a; = ujon and the troughs at a; = a;o(n + 1/2), for n = 
(a)l, (6)2, (c)3 respectively. Both the full calculated value and 
the r, A ^ a;o approximation are plotted. Parameters: ujo = 
10, r= l,A = 0,e,7 = 1/16, 



We can now use this amplitude to calculate the behav- 
ior of the SSET. If the back-action of the SSET on the 
resonator is weak, then the oscillations of the resonator 
at the driving frequency will be relatively unaffected and 
hence the drive the SSET field feels will be unchanged. 
However, the SSET does not simply respond at the driv- 
ing frequency, but has a systematic response at multiples 
of cj£), and will also act as an additional source of noise. 

Assuming that the response of the SSET at the driving 
frequency is unaffected by the back-action, the Fourier 



transform of a is given by. 

a{ujF) = 



iAoSiajp + ujd) + i^o^cTcci^F) 



i{LOF +U!q) + ^ 



i{uJF +0^0) + ^ 



(30) 



where the first fine describes the systematic motion of 
the cavity and the second gives the noise. We now see 
that our assumption that the change in the osciUation 
amphtude (i.e. the change of a{ujD)) is neghgiblc will 
be justified whenever Ad » tJof^o'fc(^^£')- ^t other 
frequencies the SSET may have a significant effect on the 
resonator. The systematic response of the resonator at 
frequency lo is just the charge response at that frequency 
multiplied by a prefactor of magnitude cl;o/((w — wq)^ + 

72 )2. 
I ex J 

We can now write the cavity noise as a function of the 
charge noise, 



Sa\.a{^) 



n^e. + ^{S!A^) + SlM) 



4:Xq 



(W - CJo)2 + 



(31) 



with the systematic and noise-induced parts of the charge 
spectrum calculated driving amplitude A. We see that 
the charge noise spectrum appears directly in the expres- 
sion for the cavity noise, and so the sidebands discussed 
in Sec. IIIII should also be present. Although the size of 
the peaks is suppressed, we find that (for n — 0), the 
ratio of the on- and off-peak noise is unchanged. 



A. Back-action damping and temperature 

In the preceding section, we assumed that the back- 
action was weak enough that the SSET damping could 
be neglected. However, it is well known that the effect of 
the back action on the dynamics of the resonator can 
be significant. In particular, as well as providing an 
additional source of noise for the resonator, the SSET 
can also act to provide an additional source of damp- 
ing [13, m, [M [lE [II- In these situations, this effec- 
tive damping will have a significant effect on the noise 
properties of the resonator even when the resonator is 
strongly driven. Furthermore, the SSET can cause the 
total resonator damping of the resonator to become neg- 
ative, and hence the resonator can be driven into a self- 
oscillating laser-like state even in the absence of exter- 
nal driving [ij, [T^ [ij, [I^. In this section we review 
how the systematic response of the SSET can act as 
an amplitude-dependent damping of the resonator, and 
present some simple analytic approximations before de- 
scribing how this influences the cavity noise spectrum. 
We show that the calculation of charge noise presented 
in Sec. Illll can be used to generalize the previously known 
expressions for the effective damping and temperature of 



a slow (wq ^ r) resonator to arbitrary frequency using 
a linear-response-like approach [27,] . 

The calculation of the damping proceeds as follows; the 
resonator amplitude changes only slowly, so we can aver- 
age the effect of the systematic SSET motion over a sin- 
gle resonator period. Following the calculation given in 
[15| we obtain an expression for the amplitude-dependent 
frequency shift and damping introduced by the SSET, 



ilssiz) + iSu;o{z))a ^ 
. XsUjQ (3r -I- 2jtt'D)e 



2Xq (r + icj£))2 

/ 1 



y^ J-n (z) Jl-n {z) 

n 

1 \ 



+ i{LUDn — A) -^ + i{ujDn + A) 



(32) 



In Fig. [Hl^a) we plot the damping as a function of am- 
plitude. The resonator A amplitude is then found by 
solving the self consistent equation for a in the rotating 
frame. 



= -i{cjQ - (jJD)a - 'fexa + iAd 
-{lss{z) + i5ujo{z))a. 



(33) 



With the resonator amplitude found, we can now go 
on to calculate the cavity noise. Before doing this, we 
look at some simple approximations to the damping. In 
the limit z — > 0, we only need include the n = 0, 1 terms 
in the sum and take the first order in z to get a linear 
damping. 



Tss 



xiLopAe^jVil?, 



^D 



(34) 



where we note that it is the driving frequency lod that ap- 
pears in the above expression rather than wq as the driv- 
ing field means the periodic motion is at this frequency. 
In the absence of external driving, lod is replaced with 
loq in the above expression. Equation (j34p then reduces 
to the previously known effective damping ^1^ [13, [lM| ■, 
but in the limit ej ^ F and extending to all resonator 
frequencies [31] rather than just wq "C T- In Fig- E^b), 
the linear damping is compared to a value obtained nu- 
merically from the mean-field equations for a range of 
values of lod/T. 

We also find that we can get a relatively simple expres- 
sion that is exact to 0{z^), and valid for finite driving 
in the limit w/j ^ F, and for a detuning that is not too 
large, |A| < loo- Including only the n = — 1..2 terms 
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FIG. 3: (color online). Plot of the SSET damping. Plot 
(a) shows non-linear damping, with —'yssA plotted against 
A. Plotted are the full expression for the damping, 7ss, the 
linear damping, ^ii„, and approximations for the nonlinear 
damping using one (71) or two (72) terms from Eq. (|35[l. 
with K = 0.1, A = -0.2, ej = 0.05, F = 1. Plot (b) shows the 
linear damping over a range of resonator frequencies, where 
the points are the damping calculated numerically from the 
mean field equations (k = 0.1, ej = 0.05, F = 1). 



gives, 



,^o2Jo(z)Ji(z) 'iAe^jxlMz)J2iz) 






(3r -f 2iu)D)3iTLL)D 



(35) 

This approximation to the amplitude-dependent damp- 
ing is plotted in Fig. ^h). 

The fluctuations in the charge act as an additional dif- 
fusion term for the resonator field [30| ■ When there is no 
external driving, we can insert the z —f 0,LUr ~ ujq limits 
of the damping and the charge noise to obtain a simple 
effective temperature. We find, 



(2nss + 1) 



'- + A' 



which again agrees with the previously known results [l6|, 
[l3,[3l| for this system. Although this expression has here 
been derived in the low-ej limit, other calculations [3l| 
suggest that this expression is exact for all ej. 

Interestingly, this expression for the effective tempera- 
ture is identical to the expression found when the me- 
chanical resonator is coupled to an optical cavity (or 
equivalent systems) rather than an SSET [32, J3j, |34| . 
with the quasiparticle tunneling rate F simply replaced 
with the optical cavity damping. This is a rather sur- 
prising result, as we have two very different systems, an 
harmonic oscillator and a 3-level SSET, providing the 
same effective temperature. The temperature is some- 
how insensitive to the details of the measuring device - 
in particular its finite level structure. This merits further 
investigation, and it would be interesting to see if other 
related measuring devices also lead to the same expres- 
sion. 



V. SENSITIVE DEPENDENCE OF THE 
SIDEBANDS ON THE DRIVING AMPLITUDE 

In this section we discuss in detail the rather sensitive 
dependence of the height of the sideband peaks on the 
driving amplitude. Due to the Bessel functions, the sys- 
tematic and noisy response of the charge at the sidebands 
oscillate rapidly as a function of z. We investigate this 
and suggest that it could be used to detect the presence 
of small changes in the amplitude of a large driving force. 

The power of the systematic peaks in the frequency 
spectrum is smaller than that in the noisy peaks by a 
factor of ej, and decreases as l/z^, so we concentrate on 
the peaks in the spectrum due to the fluctuating terms. 

We envision driving the system at some (possibly large) 
value of z = zq. As the heights of the peaks at the 
sidebands oscillates rapidly as a function oi z = zq + 5z, 
a small change of 6z will lead to a large change in the 
height of the sideband peaks. 

The sensitivity of the detector depends on the height of 
the peak at zq compared to the minimum height that the 
peak can have, i.e. the floor in the frequency spectrum 
at this point. Equation ([25|) gives an expression for the 
heights of the peaks as a function of z, and thus we have 
a simple expression for the ratio. 



R 



^Jn \^n 



2Jq{z„ 



4(7 yj [Zn 



2Jo(z„ 



(37) 
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(36) 



We find this gives a reasonable approximation for odd 
n, and for low values of z when n is even. However, 
the approximation breaks down when z becomes large. 
When z is large, the Bessel functions tend to their asymp- 
totic limit Jn{z) -^ -^ cos(z — n^ ~ j)i and we find that 
J2 — Jg for all even n, so Eq. |37|) diverges at the point 
where Jo{zmin) = 0. In this limit, in order to obtain the 
correct value for the ratio, we need to include terms to 
the next order in 1/lu]j. 



However, examining Eqs. ([2T|) and ([22)) . we see that all 
the terms in the sum are of order l/uij^, i.e. to go to next 
order requires us to perform an infinite sum. Fortunately, 
this can be done in the large z limit when the Bessel 
functions take their asymptotic form. 

For example, we can rewrite Eq. (j22|) at cj w ujone, 
where ng is even and rig > 0, separating out the odd and 
even terms to obtain, in the limit F, A ^ wo. 
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We now use the fact that in the large z hmit, Jn[z)'^ has 
the same value for all even and all odd n. This allows 
us to perform the sums exactly, using expressions like 

y^, -K — ^r '^o obtain. 



qno 



'M 



1 



12^2, 



(wDne)2(27r) \AujI, 

Performing similar calculations for Eqs. (j2ip and 
gives asymptotic forms for the charge noise for at the 
odd, even and zeroth peak to leading order. 



(39) 
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(l + sin2(z--)). (40) 



The odd peaks have a ratio between maximum and min- 
imum of order one, so will not be of as much use. The 
even peaks have maxima at z = rmr + ^, minima at 
z — rail + ^ , and are most sensitive to changes in z at 
z — '^. The ratio between the height of the peaks at 
their maxima and at their minima takes a simple form 
when n > 0. We find. 



-Rn, 



>o 



^ + A2 



(41) 



so the contrast between the peaks at their largest and 
smallest can become large in the limit uj^ ::^ F, A. 



are actually measured. In order to connect the dynamics 
of the resonator to measured quantities we need to con- 
sider how the resonator is coupled to the external world. 
What form this takes depends on what the resonator ac- 
tually is. For example, if it is the microwave field in a 
superconducting co planar cavity, we can use the input- 
output formalism of quantum optics to relate the field in 
the microwave cavity to the many mode fields in trans- 
mission lines connected to the cavity. In this situation, 
the damping and noise on the resonator is attributed to 
fields external to the cavity and these external fields are 
ultimately what is measured. In the simplest situation 
we could imagine a single side cavity with a quantum lim- 
ited input field. The output field from the cavity then 
contains a component of the reflected input field as well 
as the field transmitted from the cavity itself. This is 
the model we will adopt here as we can easily apply the 
input-output theory of quantum optics 35] as co-planar 
superconducting cavities are in the highly underdampcd 
limit appropriate for this formalism. 

Another possible realization for the resonator is a 
nanomechanical oscillator. In this case, we need an ex- 
plicit transducer model for the way in which the dis- 
placement of the nanomechanical resonator is measured. 
A typical example would be to capacitively couple the 
nanomechanical resonator to a microwave cavity [36] . In 
that case the nanomechanical resonator is coupled to 
more than one bath: the finite-temperature mechanical 
bath in addition to its irreversible coupling to the mi- 
crowave field propagating into and out of the transducer 
cavity. For a fast and efficient measurement, the mi- 
crowave cavity would be strongly damped in which case 
the resonator would see the bosonic bath due to the elec- 
tromagnetic fields on the transmission lines directly. 

The quantum Langevin equation for the field is given 
in Eq. ([5]). In the microwave realization, we will assume 
a single side cavity and that the only source of damping 
for the cavity field is in fact its coupling to the input 
and output fields at the open end of the cavity. In that 
case the noise operator, r/a is written in terms of the 
multi-mode field amplitude input to the cavity: 77a (i) — 
^lexO-init)- The output field from the cavity is related 
to the input field and the intracavity field by 



lo(t) = 



.w 



(42) 



In terms of the Fourier transformed operators, the input 
and output fields are related by 



ao{uj) = - 



i,/y^AD6{uj +ujd) 



^ - i{ujD + w) 



'-if- +i{uJD+^) 



am(w) 



(43) 



VI. CAVITY OUTPUT FIELD 

In this section we describe how the noise spectrum of 
the field in the cavity is transferred to the quantities that 



where the last term represents a phase shift between in- 
cident and reflected fleld components from a single-sided 
cavity, and 



X 



UJXa 

2^ 



(44) 
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and accii-^) is the Fourier transform of the island charge 
operator acc{t)- In general this is itself a function of the 
intracavity field and so (Jcci^^) is a complicated convolu- 
tion of a very nonlinear operator- valued function of a(uj) 
and thus Ea. (|43)) is not an explicit relation between the 
input and output field components. However if we adopt 
the approximation implicit in Eq. (|30p. we can write it 
in terms of a systematic component and noise component 



ao{uj) 






i^/l^X^^cci^) , ^ - i{^D + w) 



+ 



^ +i{LOD^ijj) ^ +i{iOD+0j) 



am(w) 
(45) 



We thus see that the noise power spectrum for the field 
output from the cavity is 



Sout{^) = / dw'{al{uj),ao{uj')) 

•J —OQ 

= n{u}) + — 3 (46) 

(w - uqY + ^ 

and we see that, for small enough n, the noise of the 
cavity is transferred to the output field, and thus the 
emitted field can be used to detect the noise in the cavity. 



VII. CONCLUSION 

We have used a Langevin equation approach to investi- 
gate the frequency spectrum of a superconducting single 
electron transistor coupled to a resonator under periodic 
driving of the resonator. The fluctuating noise terms al- 
low us to describe the correlations in the SSET due to 
the finite level structure. This approach allows the cal- 
culation of the spectrum of the charge noise in the SSET 
and the resulting effect on the resonator field in the limit 
of low Josephson energy, ej <C F. We found that the 
charge noise consists of a series of peaks at multiples of 
the driving frequency, and calculated the heights at and 
between these noise peaks in the limit of fast resonator 
oscillation. 

We have calculated the effect of the SSET on the cav- 
ity, and in particular calculated an effective amplitude- 
dependent damping and temperature that is valid for all 
resonator frequencies. We have shown that the peaks at 
the sidebands depend rather sensitively on the driving 
amplitude and show how this could provide a measure of 
small changes in the driving amplitude. 

We thank Andrew Armour and Thomas Harvey for 
helpful discussions. D.A. Rodrigues acknowledges fund- 
ing under EPSRC grant EP/D066417/1. 



APPENDIX A: DERIVATION OF THE 
LANGEVIN EQUATIONS 

In this appendix we describe in more detail our use 
of the Langevin equations and show how they can be 
derived from the master equation. We show that this ap- 
proach reproduces the results obtained by other methods 
for an SSET uncoupled to a resonator. 

We assume that the noise in the system is entirely de- 
termined by the second order correlation functions for the 
dynamical variables. We find that we can use Langevin 
equations as a useful tool for keeping track of the dynam- 
ics of the first and second moments. If the Langevin equa- 
tions give the same equations of motion for the first two 
moments as a master equation or Focker-Planck equa- 
tion then, as far as a noise calculation is concerned, the 
two are equivalent. This equivalence is discussed in Ref. 
[35| , and a description of a noise calculation given for the 
case when the conservative terms and the diffusion in the 
Langevin equation are constant in time. Here we focus 
on the situation where the external driving means that 
these terms are time-dependent. 

Although Langevin equations can be derived from a 
microscopic picture, in which the noise terms •q{t) repre- 
sent the correlation functions of the external bath (in this 
case the microscopic electron energy levels in the leads) 
here we take a functional approach and consider them 
simply as a tool to allow us to describe the correlation 
functions of the system. In this approach we find that 
we can derive the properties of the noise operators from 
the master equation. A general set of Langevin equations 
x{t) with systematic evolution described by the matrix 
A{t) and fluctuations E{t) gives, 



?(i) = ~A{t)x + E{t) 



(Al) 



x{t) = x(0)e-^(*) + e~^(*) / e^(*')^(i')clt' (A2) 



where lit) — ^ A{t')dt. From Eq. (jAip we can calculate 



equation of motion for the variance x(0 — x{t)x^{t), 

m = x{t){-x'^{t)A^{t)+E^{t)) 

+ {-A{t)x{t)+E{t))x'^{t) (A3) 



Inserting Eq. (|A2p and taking the ensemble average 
gives, 

(X(i)) = -iA{t)(xit)) + {xit))A^it)) 

+ (E{t) I /^^(i')e^"(*'Mi'|e-^"W 



E{t')dt' E^{t) ) (A4) 



e '^« / e^(*'; 
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If the dynamics of the system is Markovian (as in our 
master equation), the fluctuating noise correlators will 
be (5-correlated. Indeed, a microscopic derivation shows 
that (5-correlated operators are obtained in exactly the 
limit that the master equation becomes Markovian, i.e. 
when the decay of the correlation functions of the leads 
is much more rapid than any timescale in the system. 

Writing the correlation matrix of the noise terms as 
(^E{t)E'^{t')) — 5{t — t')G{t), we obtain an expression 
that relates the rate of change of the variance matrix to 
the Langevin correlators, 

(XW> = -{A{t){x{t)) + {x(t))A^{t)) + G{t){kh) 

which reduces to the usual [S^ ^(x) + (x)^"^ = G for the 
case of a time-independent A. The fluctuations represent 
deviations from the mean field so we can also calculate G 
by comparing the true evolution of the second moments 
with the mean-field-only evolution, which is sometimes 
more convenient in practice. Defining ^ = —A.x, the 
elements of G are also given by. 



with all other correlators equal to zero. We note that 
the commutation relations are preserved, e.g. Go2,22 7^ 
G22,02, retaining the quantum nature of the problem. 
It is also worth re-emphasizing that this functional ap- 
proach where we derive the correlators through the equa- 
tions of motion for the variance means that we are only 
capturing 2nd-order correlations. This is adequate for 
our purposes, as the noise only requires these terms, 
but these Langevin equations give no information about 
higher-order correlations. 

We can check the validity of this approach to deriving 
the Langevin correlators in a simple case where the noise 
can be calculated by other methods. For the case of an 
undriven, uncoupled SSET, inserting the Fourier trans- 
forms of Eqs. (Ill) to ([9]) along with the expressions for 
the correlators from Eq. (|A7|) in the expression for the 
noise, Eq. ((M]), gives an expression that is mathemat- 
ically identical to the noise as calculated more directly 
using the expression 



G 



ij — Y^i-^j) 



k^iX 



l-^J I 



{X^^j) 



(A6) 



Using either Eq. lASl or Eq. IA6| we find that the correla- 
tors for the SSET are given by, 

Goo, 00 = —Goo, 11 = —Gil, 00 = r(o'ii(t)) 

G22,22 ~ —Gil, 22 = — G22,ll = r(cr22(i)) 
Go2,22 ~ — Go2,ll = r(cro2(i)) 
G22,20 ~ —Gil, 20 = r(cr2o(^)) 

Gii,ii = r((aii(i)) + (a22(i))) 

Go2,20 = r((aoo(i)) + (aii(i))), (A7) 



(cr,j (t + T)aki (t)) = (exp(-ylT)cr„- {t)aki (t)) (A8) 



which can be evaluated by exactly diagonalizing A. In 
particular, this is true for all values of F compared to 
the other SSET timescales. A, ej. We can also use the 
converse of this argument to confirm that the noise terms 
must indeed be ^—correlated. 



[1] Y. Nakamura, Yu. A. Pashkin and J. S. Tsai, Nature 398 [11] 

786 (1999) 
[2] J. M. Martinis, S. Nam, J. Aumentado and C. Urbina, [12] 

Phys. Rev. Lett. 89 117901 (2002) [13] 

[3] A. Cottet, D. Vion, P. Joyez, P. Aassime, D. Esteve and 

M. H. Devoret, Physica C 367 197 (2002) [14] 

[4] J. E. Mooij, T. P. Orlando, L. Levitov, L. Tian, C. H. 

van der Wal, and S. Lloyd, Science 285 1869 (1999) [15] 

[5] A. Blais, R.-S. Huang, A. Wallraff, S. M. Girvin, and R. 

J. Schoelkopf, Phys. Rev. A 69 062320 (2004); A. Wall- [16] 

raff, D. I. Schuster, A. Blais, L. Frunzio, R.-S. Huang, 

J. Majer, S. Kumar, S. M. Girvin and R. J. Schoelkopf [17] 

Nature 431 162 (2004) [18] 

[6] M. P. Blencowe, Phys. Rep. 395 159 (2004) 
[7] D. Mozyrsky and L Martin Phys. Rev. Lett. 89 018301 [19] 

(2002) 
[8] A. A. Clerk and S. M. Girvin Phys. Rev. B 70 121303(R) [20] 

(2004) 
[9] J. Wabnig, D. V. Khomitsky, J. Rammer and A. L. She- 

lankov Phys. Rev. B 72 165347 (2005) [21] 

[10] D. Mozyrsky, L Martin and M. B. Hastings Phys. Rev. 

Lett. 92 018303 (2004) [22] 



A. D. Armour, M. P. Blencowe and Y. Zhang Phys. Rev. 
B 69 125313 (2004) 

R. S. Knobel, and A. N. Cleland Nature 424 291 (2003) 
D. A. Rodrigues and A. D. Armour New. J. Phys. 7 251 
(2005). 

D. A. Rodrigues, J. Imbers and A. D. Armour Phys. Rev. 
Lett. 98 067204 (2007) 

D. A. Rodrigues, J. Imbers, T. J. Harvey and A. D. Ar- 
mour New. J. Phys. 9 84 (2007). 

M. P. Blencowe, J. Imbers and A. D. Armour New J. 
Phys. 7 236 (2005) 

A. A. Clerk and S. D. Bennett New J. Phys. 7 238 (2005) 
S. D. Bennett and A. A. Clerk Phys. Rev. B 74, 
201301 (R) (2006) 

M. D. LaHaye, O. Buu, B. Camarota and K. C. Schwab 
Science 304 74 (2004) 

A. Naik, O. Buu, M. D. LaHaye, A. D. Armour, A. A. 
Clerk, M. P. Blencowe and K. C. Schwab, Nature 443, 
193 (2006) 

P. Filipowicz, J. Javanainen and P. Meystre, Phys. Rev. 
A 34, 3077 (1986). 
H. Walther, B. T. H. Varcoe, B.-G. Englert, and T. 



12 



Becker, Rep. Prog. Phys. 69 1325 (2006) 
[23] O. Astafiev, K. Inomata, A. O. Niskanen, T. Yamamoto, 

Yu. A. Pashkin, Y. Nakamura, and J. S. Tsai, Nature [31] 

449 588 (2007). [32] 

[24] M.-S. Choi, R. Fazio, J. Siewert, and C. Bruder, Euro- 

phys. Lett. 53 251 (2001); M.-S. Choi, F. Plastina and R. [33] 

Fazio Phys. Rev. B 67 045105 (2003) 
[25] A. A. Clerk, S. M. Girvin, A. K. Nguyen, and A. D. Stone [34] 

Phys. Rev. Lett. 89 176804 (2002) 
[26] J. Hauss, A. Fedorov, C. Hutter, A. Shnirman and G. [35] 

Schon Phys. Rev. Lett. 100 037003 (2008) 
[27] A. A. Clerk Phys. Rev. B TO 245306 (2004) [36] 

[28] F. Marquardt, J. G. E. Harris and S. M. Girvin Phys. 

Rev. Lett. 96 103901 (2006) [37] 

[29] M. Ludwig, B. Kubala and F. Marquardt 

\arXiv: 0803.3714 (2008) 
[30] Ya. M. Blanter, O. Usmani and Yu. V. Nazarov, Phys. 



Rev. Lett. 93 136802 (2004); ibid., Phys. Rev. Lett, (er- 
ratum) 94 0409904 (2005) 

A. A. Clerk, unpublished; D. A. Rodrigues, unpublished 
I. Wilson-Rae, N. Nooshi, W. Zwerger and T. J. Kippen- 
berg Phys. Rev. Lett. 99 093901 (2007) 
F. Marquardt, J. P. Chen, A. A. Clerk and S. M. Girvin 
Phys. Rev. Lett. 99 093902 (2007) 

M. P. Blencowe and E. Buks Phys. Rev. B 76 014511 
(2007) 

D. F. Walls and G. J. Milburn, Quantum Optics (Berlin: 
Springer- Verlag) (1994) 

M. J. WooUey, A. C. Doherty, G. J. Milburn, K. C. 
Schwab, arXiv:0803.1757 (2008). 

C. B. Doiron, B. Trauzettel and C. Bruder, Phys. Rev. B 
76 195312 (2007) 



